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CHAPTER 1 


INTRODUCTION 

1.1 Application of Stochastic models in fertility 
has been an important component of investigation in 
demography concerning the human reproductive process. 
Some researchers have shown their interest to develop 
probability models for the number of conceptions in a 
certain period giving rise to discrete probability 
models (Singh 1963), (Pathak 1966). A good amount of 
research has been carried out related to probability 
models of the waiting time distribution for the first 
conception as well as between conceptions of different 
orders . 

The researchers involved in the field of 
demography have utilized opportunity to develop 
analytic model building utilizing the theory of 
Stochastic Processes and Renewal Theory. Family 
building models are particularly concerned v.'ith the 
series of conceptions, births etc. that constitute the 
reproductive history of a couple. 
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Another aspect of this type of research in 
this area is to find practical solutions to problems 
connected with decision making and evaluation of family 
planning programmes. The problem of measuring the 
impact of family planning programme on fertility has 
become a major research issue for demographers as well 
as policy makers. Mathematical demographers have 
contributed considerable literature regarding fertility 
analysis with Stochastic models. These models have 
advantage over the usual deterministic models. For 
example, with stochastic models in fertility analysis 
one can obtain sample estimates of the demographic 
parameters with exact specification of Sampling errors. 
Considerable researchers have worked out the 
probability of conception allowing for the non-fecund 
period of pregnancy. Models have generated to estimate 
these probabilities (some -researchers have taking 
foetal wastage into consideration on the other hand 
some have ignored by assuming complete pregnancy 
leading to live birth). These research v.'ork are based 
on the data source providing the distribution of birth 
intervals betv/een tv/o successive birth orders. Mitra 
and [banner jee (1982) v.liile estimating the probability 
distribution between tw'o successive live births also 
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investigated the same distributions in a different way 
to estimate the lengths of the non susceptible period 
that begin from the time to conception. 

' Since conception can occur during only a 
small segment of the ovrian cycle and mean length of 
the cycle is about 30 days thus the conception time was 
treated as i) Discrete ( each unit being taken equal 
to one month) 

ii) Continuous. 

The variables under consideration in such 
probability models are thus the random lengths of 
intervals between events. The family building models 
were thus introduced in both discrete and continuous 
time considering both hetrogeneous and homogeneous 
population. The importance of stochastic models in 
fertility analysis has been encountered by Talwar 
(1966) . By considering t!ie distribution of the time 
interval between consecutive conceptions he has 
obtained the joint distribution of the period of post 
partum amenorrhea, gestation period and fecundable 
period. From this model one can estimate the parameters 
concerning the period ot post-partum araenorrhoea 
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(P.P.A.) from the data of interlive birth intervals 
which are easily available from most of the 
retrospective demography surveys. 

Biswas (1975) has shovjn that with application 
of stochastic models in fertility analysis one can 
estimate proportion of women exposed to the risk of 
fecundity, the proportion under family planning 
practice and from which the proportion of births 
interrupted due to family planning practice can also be 
estimated. Stochastic models can also be used in 
analysing demographic data based on follow up surveys. 
In such surveys either the period within v/hich data is 
collected is taken as fixed so that the proportion of 
observations under consideration during that period is 
a discrete or continuous random variable within which a 
fixed number of observations under consideration lie. 
There had been a study with a ne’w investigation which 
nuiy highliqlit the interrelationsiiip between two useful 
variables in the human reproduction process; viz the 
period of post partum amenorrhoea and 'the lactation 
period which are useful items of consideration for 
averting methods which has been done by the application 
of a bivariate survival model due to hreund (1961). 
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1.2 REVIEW OF EARLIER STOCHASTIC FERTILITY MODELS: 


The application of stochastic models in human 
fertility started in the year 1924 with C. Gini's paper 
presented at the International Mathematical Congress in 
Toronto. In recent times considerable interest is 
generated in the analysis of pregnancies of cohorts of 
married women. Henry (1953) has emphasised on the 
patterns of live births. The human reproductive process 
is thus compared with Geiger Muller Counter where every 


event occurring with 

an 

intensity 

is follo\i/ed by a 

dead time (or blocked 

time) 

during 

v;hich 

no further 

success is possible. 

G. 

M. 

counter is 

classified 


according as the dead time is fixed or random variable. 
The term "counter is paralyzed" is used in this respect 
which is comparable to the non - susceptible period 
following a pregnancy in the human reproductive 
process. The type II G. M. counter are contrasted from 
type I in the sense that not only ttiere cannot be any 
success during the dead time in type II during the dead 
time further prolongs the dead time without registering 
the event during tlie blocked period. 




For the analysis of human reproductive 
process, counter model I with fixed or random variable 
dead time appears to be more appropriate. Human 

reproductive process are directly derivable from the 
existing results of G. M. counter models. Further the 
renewal theory technique is amenable for the 
development of probability models on fertility and the 
human reproductive process. The models can be 
classified in two categories; 

1. Deals with the waiting time distribution between two 
consecutive events. 

2. Deals with the probability distribution of the 
number of events during a fixed time. 

In first category time is a random variable 
and the system conforms to ordinary or delayed Renewal 
process (Non - Markovian in character) . In the Second 
category the number of events is a random variable 
during fixed time and conforms to a renev/al counting 
process. 

In India, the making of Stochastic model 
representing the human reproductive process V\?as first 
done by Dandekar (1955). Where Dandekar assumed one to 



one correspondence betv/een pregnancy and time birth and 
a constant time independent Poisson intensity for 
births/ pregnancies of all others. A fixed non- 
susceptible period ( dead time) was assumed following 
every pregnancy/ birth in the pattern of counter Model 
type I. The probability distribution of the number of 
events in fixed marital duration (o,t] has been 
obtained from the Modified Poisson Distribution as 
evolved by Dandekar. The distribution considered is 
Poisson Renewal Process where every renewal is followed 
by a dead time of fixed measure. Thus it is obtainable 
as an example of Counter Model Type I with fixed dead 
time as considered by Takacs (1960), Dharmadhikari 
(1964), Biswas (1973). Basu (1955) adopted difference 
equation techniques which provides alternative 
derivation of asymptomatic pregnancy rate. 

Brass (1958) obtained Negative model from 
Dandekar's distribution ignoring tlie non-susceptible 
exposures follov.dng every birth and obtained a 
theoretical model for the probability distribution of 
births. However, the model did not give good fit in 
empirical data as infecundable exposure following every 
birtli v.'as ignored, the probability of zero number of 
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births to mothers v/ith completed reproductive span as 
obtainable from William Brass Model did not show good 
correspondence with that of the empirical frequency for 
the same class because of presence of considerable 
proportion of Biologically sterile mothers ( or couples) 
as well as a good proportion of mothers (or couples) 
who wanted to remain effectively childless without 
being biologically sterile. 

Singh (1963) independently obtained two 
different interrupted probability models on similar 
assumption v;ith the problem of estimating the number of 
schools of fishes in a fishing ground. For example, the 
dead time "h" following each conception was comparable 
to the time for which further counting of fishes was 
temporarily abandoned following the sight of a school 
of fish and exactly "h " hours were spent for fishing 
the school. In both the models Singh ■’ assumed the 


presence 

of (i- \) of 

the couples 

to be 

effectively 

ster i le . 

V\’hile in the 

first model 

Singh 

assumed the 

montli ly 

probability of 

concept ion 

to be 

constant, in 

the second model he 

assumed the 

same 

admitting a 


probability distribution of B typt; following Potter and 


Parker . 


Another aspect of S. N. Singh's work in this 
direction comprises of analysing the building up 
appropriate probability model for explaining the random 
changes causing variability in the- number of 
conceptions during a fixed fecundable period. In the 
first course of his study Singh (1963) devoted to this 
analysis on complete conception only and obtained BAN 
estimates of the parameters of the model, whereas in a 
generalized model built up later, Singh and Bhatacharya 
(1970) admitted incomplete conceptions too in the model 
which led to several more parameters of the model. A 
partial estimation of some of the parameters was made 
while assuming the knowledge relating to the rest. 

Dharmadhikar i (1964) while attempting to 
generalize Dandekar's model , first of all developed a 
model based on Mixture of Exponential distribution with 
the origin shifted by K ( the length of the non- 
susceptible exposure ). Besides this generalization of 
Dandekar's model modified into a mixture of _ two 
distributions, Dharmadhikar i developed a generalized 
time related dependent process (Doubly Stochastic 
Poisson Process) and considered the application of the 
same in further models. 
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S. N. Singh (1964) further obtained the 
interrupted probability models for couple fertility 
under the assumption of the presence of a fixed 
proportion of effectively sterile couple as well as 
fixed non- susceptible exposure following every 
conception and obtained BAN estimators of the 
parameters of the model by using Central Limit Theorem. 
Modified Minimum OG ^ method has been applied to 
estimate BAN estimators. As noted by Sheps and Menkin 
(1969), Singh's model is a modified restatement of 
Dandekar's model using a Neyman's Method (1949). 

Perrin and Sheps (1964) made an explicit use 
of renewal Theory for building of probability models on 
human reproductive process that allow a variety of 
pregnancy outcomes (live birth, still birth, foetal 
wastage, spontaneous or induced abortions) and general 
distribution for durations of pregnancy and the post- 
partum non- susceptible status. The states of the 
reproductive process in Sheps and Perrin consist of : 

Sq = Non - pregnant fecundable state 


= Pregnant state 


5 2 = Post -parturo infecundable state associated with 

abortion or early foetal loss 

5 3 = Post -partum infecundable state associated with 

still birth 

= Post -partum infecundable state associated with 
with live birth 

The length of stay in each state is viewed as 
a random variable, the waiting between two renewals of 
the same state and the number of renewals of each state 
during a fixed term (o,t3 conforms to a Semi Markov 
pattern. To study the fertility process it is important 
to know number of passages into a specified state , the 
fertility rate (or the pregnancy rate per unit of 
time) ; and the probability distribution of the renewal 
of different states of the process in a given period of 
time. The results obtained from these models have 
implication both for the analysis of data on the 
distribution of births to a group of v;omen in a 
calender period and for the planning and evaluation of 
family planning programmes. In this direction a family 
building model has been set up by Sheps , Menkin and 
Radick (1969) . A very important model for studying the 
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short term effects of change in reproductive behaviour 
( Sheps and Menkin (1971) ) on intrinsic birth rate and 
a model for concept ive delays in a heterogeneous group 
have also come up. 

Theoretical investigation for developing 
Inter-arrival time distribution under weighted Poisson 
Process by applying Palm Probability technique has been 
done. The renewal intervals corresponding to different 
order of conceptions become correlated because of 
variability between individuals thereby distributing 
the entire renewal structure. Multistage Markovian 
chain model with forward and backward transition 
probabilities represented by a density dependent birth 
and death process has been developed to study 
effectively the efficacy of a particular sterilization 


policy 

motivated to 

reduce the birth 

rate in 

a 

population. ■ 

Attempts 

have been 

made to 

. analyze 

the 

direct 

and 

indirect 

strategies 

for the 

reduction 

of 


conception rate. While classifying family limitations 
methods as direct and indirect , the inter-live birth 
intervals bet\>/een 1-2, 2-3, & 3-4 have been vjorked out. 
G. M. counter technique has been applied to take 
account of the infecundable period. There has been a 
nev; application of survival analysis technique for 


obtaining inter live birth intervals by employing a 
general hazard model which can take account of any type 
of fertility behaviour which one may theoretically 
assume ( Pachal 1994) . 

1.3 A SURVEY OF SOME OF THE METHODOLOGIES EMPLOYED; 

1. Palm probability is defined as the conditional 
probability of a specified number of events in a time 
interval given that an event has happened at the 
beginning of the interval. Cox and Isham (1980) 
describe Palm Probability as follows : 

For u > v, let N(u,v) be a random variable giving the 
number of events occurring in (u,v) and be the 

sequence of the intervals between i th successive 
events (i= 1,2,...) in a process starting from an 

arbitrary point. 

Consider the survival function 

K (x) = P(X>x), where X is a r.v. representing the 
waiting time for first renewal 

lim 

= 6>o+ P f N(0,x) - 0| N (-6-,o) > 0] (1.01) 
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( subject to orderliness which implies not more than 
one renewal in the infinitesimal interval ) which is 
the limiting probability that, given an event occurs 
immediately before the origin, the next event of the 
process occurs after that instant ( i.e the system will 
survive more than a period x ) . 

Now, by stationarity condition 

P[N(0,x) = O n N(-(S,0) > 0)] 

= P[N(0,x) = 03 - P [N(-5,x)=0] 

= P[N(x) = 03 - P [N(x +5)=0] (1.02) 

Now P[N(0,x) = 0 I N(-5,0) > 0)] 

P(N(0,x) =0 N(-<S,0) > 0) ] 

P[N(-<S,0) > 0) ] 

P[N(x) = O] - P [N(x+5) = 0] 

P{N(6') > 0) ] 

P [N(0,x) = 0 j N ( -5,0)>0] S~^ P[N((5>0] 


= 5 ^ P[N(x )=03 - P[ M(x+5)=0]] 

= - [P[N(x+5)=0] - P[N(x)=0]] (1.02') 

Define 

lim cS"! P [N(5)>0] = A (1-03) 

5-0 
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as the occurrence parameter A 
assume as finite. 

Further, denoting 


of tlie process v/hich i 


P[N(X) =k] 


Pk = , 1 


as the distribution of N(x), 
(X) = P(X>x) 


(1-04) 


then in the limit as (5-»o + 


= 1 im 


P[N(0,X) - 0| N (-6', 0) 


> 0 ] 


(1.04 ' 


NOW P(N(0,x) =0|N(-i,0) ,0, 


0) 


P[N((S) >0 ] from (1.04' 
Since the process is orderly , if k>0 
Pj^(x+<S) = p 


as 5-»o + 


as (S-*0+ 


o that 


- PfN(-(.V, 0) =0, 




h(u,x)=k-l] + 0(6) 


= Pk(>^) 
+ 


P[N(-^,0)>0, N(0,x)=k] 

I [ N ( - cS , 0 ) >0 , N ( u , >; ) =k - 1 ] +0(5) 


Pj,. (x+ 6 ') - pj^. (X) = 


P[ N (-5,0) .-0 , N (0 , X) =k 1 



= - P(N(-i, 0 )> 0 . P[N( 0 ,x) = k|N(-S,0)>0] + 

P[N(-i,0)>0. P[N( 0 ,x).k-l|N(-s,o,> 0 ] + 0(J) 
Thus ^-1 [P(N(-^,o,>o). P[N( 0 ,x,=ki N(-i,0,>0) 

- - PtN(-J,0)>0). P(N(0,x).k-llN(-i,0)P01 + 0(1) 

Xx (X) A - Dj, (P^ (X)] as i -0 .... ( 1 . 05 , 

where Dj, denotes the derivative. 

The equation 1.05 links the distribution of 
the interval between suocesssive events with survivor 
tunotion K, ,x, to that of forward recurrence tine with 
survivor function P^fx, (where P„(x, is the probability 
that starting fron an arbitrary tine instant, there are 
no events in the following interval of length x. 

1-4 PALM'S INTEGRAL EQOATION FOR PALM PROBABILITY AS A 

limiting probability : 


qcnctal iGsults connecting distributions 
ot events conditional on a point at the origin with 
those Where the origin is an arbitrary instant nay be 
Obtained. Assuning that the process is conpletely 

ionaiy and has a tinite occurrence parameter y and 

is orderly so that y is enmi +- g-, 

y IS equal to the rate of the 

process. Then ior eacl, >; > o the - i i- 

/ {jcilm distribution 


is a discrete distribution defined 


M = lin^P , , I ^ 


for k= 0,1,2 


in a careful mathematical development the 

existence of v n,nH „ / > 

^ ^ more generally 

of other limiting nroh^hi' i i -i- • 

g probabilities Of events B given 

N (-<S.0)>0 of the form 


lim P [ B i N (-^,0) > 0 ] 


5 -► 0 + 


( 1 . 07 ) 


has to be proved. 


Let the probability measure II be defined for 

events B on those processes which have , ■ . 

wjiicn have a point at the 

origin. The measure II can then be shown to satisfy 
" P t B I N (-.,0) > 01 


<S - 0 + 


for a wide class of events r m, 

s B. The measure II is 

called the Palm measure of the process me 

fiuoess. In equation 

(l-oa, the distribution of the interval measured from' 
an arbitrary time instant to the newt point of the 


process is linked to that of t-i 

^ interval between 

successive points. In i-he 

the same v/ay tlie functions 
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77]^ (x) defining the palm distributions given in (1.06), 
which specify the distribution of the number of events 
in an interval of length x given a point at the 
origin, can be connected with the functions Pj, (x) 
which give the distribution of the number of events in 
an interval of length x starts at an arbitrary 
origin. These connecting equations are known as 

Palm-Khintchine equations and may be derived as follows 

[Pjc(x+<5) - Pk(x)] = -S~^ [P[N(-5,0)>0]]. 

[P(N(0,x) )=k| N(-<S,0)>0)] 

hence for limit as £-*0+ and for right hand derivative 

Dx 

[Pk(>^)] = P[n(-5,0)>0] . [7rj^_^(x)] 

= (from (1.03)& (1.06) ...(1.09) 

The corresponding equcition for k = 0 has 

already been derived and is 

DxfPo(^^ ^ " “A (1-10) 

It follows from (1.09) and (1.10) 

-1/ X Dj.. [Po(x) +. . . +pj,, (X) ] - 
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-1/ A D^[P[N(x) < k]] = 7rj^(x) 


Therefore, the probability of having exactly 
k events in (0,x) starting from an event at 0, can be 
obtained by differentiating the probability of getting 
not more than k events in (0,x) where 0 is an random 
time instant. Alternatively from (1.11) and 1.12). 

Po(^) +Pi(5<) + + P}^(x) = P(N(x) < k) 

= 1- A P ^ K (1.13) 


So that the probability of getting not more 
than k events in (0,x), where 0 is an random instant 
can be obtained by integrating the probability of 
exactly K events in the interval when there is an 
event at 0 . In addition the right hand side of (1.13) 
is equal to 


1 " A / "k du = A J (u) 


du 


Pn(^) +Pi (>-) +• 


■ • - + Pk ( 1 ~ Aj k ( ’'A du 

0 

OO 

Pq (‘^’) + Px(“o) + +Pj,. (uo) = 1- ^ J du 

O 

^OO 

--> A du = 1 


o 


■7C 


1 - A 


r 


ta 




du = 


A y 77j^(u) du - a 

£j o' 


-n-j^(u) du 


= A / TTj^Cu) du 


OO Cfi 

Again Aj du = n-^ y+x dy putting u=y+x 

0 

Hence Po(>^) +Pi(x) + + Pj^Cx) = 

P|N(x) < k] = a/’^ (y+x) dy (1.14) 

0 

(1.14) may be justified by the argument that 
if 'O' is an arbitrary time instant and there are no 
more than k events in {0,x) then there most exist an 
event with coordinate y, for some Y > 0, such that 
there are exactly k events in (-Y, x) . Since the 
process is orderly, the probability of an event in 
i~y I "Y ti^) is A + 0(<5) and therefore 


P[ N(x) < k] 


j '"k 


(y+x) dy 


The equations (1.09), (1.10) and (1.13) can 
be summarized by using probability generating 
functions . 
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For if it is defined 

G(Z,x) = S (X) 

Go(Z,X) = f Z^ 77j, (X) 

l<Z.o 

so that G refers to an arbitrary origin while Gq refers 
to the situation given a point at origin. 

In fact, 

(X) = G (X) 

C~K. 

where the right hand side denotes the cumulative 

distribution of [ + ^ 1 ^- 3 

obtained by k - fold convolution. 

That is 

TTj^ (x) = X - G (x) 

and the P}^(x) are given by 

■3C 

Pk(x) = - J [TTj^Cu) - du for k>0 ....(1.5) 

® 7t 

Pq(x) = 1 - (1*6) 

0 

1.5 BIRTH AND DEATH PROCESS : 

One of the obvious generalization of the pure 
birth processes is to permit X (t) to decrease as well 
as increase, for example by the death of members. These 
are more relevant to biological population in which 
both births and deaths occur. The birth and death 
process developed here is applied to some problems of 
policy making in sterilization in chapters 5 of the 
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thesis . 


Let X (t) = Size of the population at time t 
with the initial size X(0) = Kq 

and P (t) =0 where X (t) = k 

Given X (t) = k v.'e assume that the probability 

of exactly one birth in the interval (t , t+<St) is 
Ao(t) + 0 ( 6 ) 

(i) The probability of exactly one death is 

M,^(t) + 0(5) 

(ii) The probability of more than one change is 0(5) . 

Therefore, the probability of no change in 
(t, t+5t) is 

1 - A)<(t) - Mk(t) +0(5) . 

Consequently, the probability Pj^(t + 5) at 

some time t+5 may be given by the kolmogorov equation 

Pj,,(t + 5) = pj,(t) [ 1 - (AK(t) +/.4j,(t) +0(5))] 

+ P).;_i(t). Ak-i(t).5+ pj.^^(t). .5 

+ 0(5) (1.17) 

This gave us the corresponding system of differential 
equation 


dPo(t) = - [Ao(t)+Mo(t) ] Po(t) +Mi(t)P3L(t) .. (1.18) - 

so on till 

dPK(t) = -[Ak(t)+Mk(t) ] p^(t) .. (1.19) 

for k>l 

Under initial condition 

Pj.o(0) = 1 and P]^(0) = 0 for k=kQ 
(1.19) completely determines the probability 
distribution P}r(t) . 

Appropriate assumption may be made regarding the 
function A]<(t) and jUj^(t) to obtain stochastic 
process corresponding to empirical phenomena. 

A particular case : linear growth 


V'*-' '» 


=> d Po(t) = MPi(t) (1.20) 

dP]^(t) = - k (k+M)P]^(t) + (k-1) XPk_l (t) 

d-t 

+ (k + l)MP}^+l(t) (1.21) 


Employing the method of probability generating 
function (p.g.f.) to solve the same 


G^(s,t) = E s p.(t) ... (1.22) 

k-z-O 

From (1.21) it follov;s 

^G^(s,t) + (l-s)(A s-M ^j^(s,t) = 0 .... (1.23) 

d-t 


With the initial condition t = 0 

.^ho 


G>,(s,0) 


(1.24) 


■ 2 3 


The auxiliary equations are 


dt = 


ds 


(1.25) 


(l-s) (^s-M) 


and 4-Gx (s,t) = 0 




For \ = ^J. ■ye may use partial fractions to rewrite 
the first auxiliary equation as 


dt 


+ " is 

(A - M) (As - M) ( A-M) (l-s) .. (1.26) 


=> 


(A"M)dt = S log 


A s - M I 


1 - s 

Integrating both sides of (1.27) => 


l-s 


s - jj, 


.(A-m) t 


= constant 


(1.27) 


( 1 . 28) 


The second auxiliary equation of (1.25) 


(s,t) = constant . (1.29) 

Therefore the general solution of (1.23) is 


G^(s, t; 


0 


l-s 


As-m 


(A-m) t 


(1.30) 


where 0 ,is an arbitrary differentiable function. 
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Using the initial condition, it can be seen 


at t = 0 


0 


1-s 


As-M 


= s' 


at least for all s with is| < 1. 


0 ( 0 ) 


Hence for all 0 such that 
ko 

1 

1+6 jj. 


1+0 


Letting 1-s 

0 = 

As - M 


( 1 . 31 ) 


1+0 < 


( 1 . 32 ) 


( 1 . 33 ) 


The particular solution for the case of iJ, 
given by 

r 

(Xs-ji) + ;u(l-s) 


Gy (s, t) 


{Xs-jx) + A(l“S) e 


Putting 


for IJ. = X 

1 - 


a(t) = jx 


0 - A s 


(A-m) t 


(^ 


A 


B(t) = — a (t) 
M 


The p.g.f. can be rev;ritten as 


I 1+0 1 


A 

ko 

( 1 . 34 ) 

. 35 ) 
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Expanding the p.g.f. 


Gjj (s : t ) 


Pk(t) 


Pk(t) 


P[ X(t) = k] can be written as 



- ’"o- 


- k^+k-j-ln 

s 

^ j ^ 


1 


J 

+ (a(t)) (B(t))*'^"^ [1-Q:(t) -B(t)]^ 


for k > 1 (1.37a) 

and Po(t) = [ a(t) (1.37b) 


The use of Martingales based on a density 
dependent birth and death process is done in chapter 5 
on sterilization policy. 

A Stochastic Process : n = 0, 1, 2, 

3,... is a martingale if for all n = 0, 1,2,3, 

i) E { 1 Xj^l } < CO 

ii) E { I Xq, X^, ... , X|.j } 

(1.38) 

hold 

A more general definition : 

Let{ Xj.j } : n = 0,1, 2, 3, and (Y : n 

= 0, 1, 2, 3, ) are stochastic processes. 

We say { X^^ } is a Martingale with respect 

2G . 


to { } if for all n = 0,- 1, 2, 3, 

i) E { 1 X^l } <00 

ii) E { I Yq, Y^, ... , ^ ~ ''^n 

(1-39) 

One may imagine ( Yq, Y^^, Yj^) as 

the information history upto the stage n. 

1.6 WALD'S MARTINGALE : 

Let Yq = 0 & 

E [ exp { A Y }^} ] existing for some \ = 0 

where (Vj^) 's are independent identically distributed 
random variables for k = 0, 1, 2, 3, ... then 

X n = ( 0(A))“" e A(V. +...+ V„) 

is a Martingale with respect to Yj^. 


1.7 SUB MARTINGALE AND SUPER MARTINGALE: 


Let { ; n = 0,1,2,...} and { Yj^ : n = 0,1,2,...} 

be stochastic processes then (Xj-^) is called a Super 
Martingale with respect to (Yj^) if for all n 
i) E { I Xj } > - 00 

where x" = inf { X,0 } 


ii) 

iii) 


Vn/ Yn ] 


Xj^ is a function of Yq,Y^ 


Similarly (X^^) is called a Sub Martingale with respect 
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to (Y,^) 

if : - 



i) 

E {X/} 

< 00 v;here X"^ = Sup (X, 

0) 

ii) 

E [ X^^p 

1 Yq. Vp, Y 2 ... , Y^ ] 


iii) 

is a 

function of •' 

■ '^n 


1.8 DEFINITION OF MARKOV TIME : 

It occurs in many Contexts. A random 

variable T is called Markov time v; . r . t . Yn associated 
with a Martingale (XrjO for every n with respect to 

marker chain (Yj^) if T takes the value 0, 1, 2, 

and if for all n = 0, 1, 2, the event T = n is 

determined by (Yq, Y^, Y^^) . 

By the term 'determined' implied that the indicator 
function of the event {T=n) cam be v/ritten as a 
function of Yg, Y^, Y 2 , , Yj^. 

1.9 OPTIONAL SAMPLING THEOREM : 

Suppose (Xj-,) is a Mcirtincpile and T is a 
Markov time v;ith respect to (Y ) tlien we establish 
E(Xo) = E ... (1.4U) 

If T < 00 then lim Xrppj.^ = X rj, ; 

.ictually X,pj.|j.j = X,j, '.-.'.here n> 'P . 


2 o 




■ Thus whenever we can justify the interchange 
of limit n -> 00 and we can formulate the following : 

E[Xq] = lira E[X.pp^] = E [ lira X^^^] = E(Xy) -• (1.41) 

W-'ioo h 

l./o OPTIONAL STOPPING THEOREM: 

Let ^ Martingale and T be a Markov 

time. If 

i) p[ T < 00 ] = 1 

ii) E[ ! X^ ! ] < 

iii) lira E[ X I(T>n) ^ " ° 

EiXrj.) = E(Xq) (1.42) 

where I ^ rp v, represents the indicator function 
for T > n. 

Counter Theory is employed to evaluate the 
direct and indirect strategies for the 

reduction of conception rate. 

A Geiqer-Mu 1 ler Counter Model I is a 
registering mechanism for detecting the presence of a 
radioactive material arriving at the counter, but 
because of inertia, the counter will not register some 
of the impulses. More specifically , suppose that an 
impulse arrives at a fixed time t the counter registers 
the impulse. Tlie registration causes a dee-id time say of 
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length ir ^ and impulses during the dead time will not be 
registered by the counter. In general, the first 
impulse to arrive after the termination of the dead 
time again will be registered by the counter and 

this again causes a dead time of length say 712 
and so on. However, in the simpler cases, the impulses 
which arrive during a dead time do not cause any dead 
time so that each dead time is caused by a registered 
impulse. This is called a Geiger-Mul ler Counter of type 
I. (Pyke 1958). In a counter of type II (Smith 1958) 
each arriving impulse causes a dead time so that 
arrivals during a period of dead time prolong further 
the dead time. A counter model although basically a 
description of certain physical processes can be 
applied to many Biological, Social and Industrial 
processes . 

1.11 REPLACEMENT MODELS: 

We illustrate the replacement model decid time 
using the follovv'ing problem : 

A conception takes place with intensity A 
subject to the condition that every conception is 
followed by infecundable exposure x (fixed dead time) 
during which no further conception takes place. Then 

1 u 


to obtain the probability distribution of the number of 
conceptions in (o,t]. 


Let ’’i '^2 '^n renewal times 

(waiting time of conceptions) v/ith 'conception rate 
(Poisson intensity) having negative exponential density 
function . 


/(t)= >»e^^O<t<>x);;>>>0 


P[T^<X] = l- e”'^^ ; A > 0 ; 0< X < 00 
P [ < X } = 1 - for X>77 

= 0 otherv;ise (1.43) 


Let W (t,n) = Probability of not more than n events 
(conceptions ) up to time t 
W(t, n) = P > t] 

= - 1 - P < t] 

l-Fn+l'(t) 

Pn+1 (t) (1-44) 

where F,-i+]_(-) is the c.d.i. of tlie random variable 
and (t) = l-Fj^_j_^(t) is tiie corresponding survival 


function . 


Denoting 


L ( . ) as the Laplace transform 


L[W(t, n)] = J (l-F^^^(t) dt 


0 

1-- 1 L [F^^^(t) ] 
s s 


(1.45) 


L{F^^,(t)) = 1 (t) ] v.'here 

'} interval density function of (n+1)^*' renewal (order of 


conceptions) . 

We have f 


F 1 * { F ^ ^ * 


n + 1 

where f^ is the density function on r^, f is the 
density function of r > 2, and * stands 

for convolution; and f(n)* s n fold convolution of f ^ 

/ 


F [Fn+i(t) ] 
Now L [ f ^ (t) ] 


L[{f^(t)} {L(f(t)}^] 

A 


A + s 


and L [ f ( t) ] 


L t^n+1 (t)] 


je-st ^ e-A(F-«)dt 

7 / 

A e 

s + A 


—ns 


A 


ni- 1 


(1.46) 


(S A ) 


n + 1 


,-n'ffA 


( 1 . 47 ) 


Putting (1,46) and (1.47) in (1.45) 


L [W (t, n) ] 


S (S + X) 


(1.48) 


By taking Inverse Laplace Transform 

t e~ A (u-n77) ^n 

W (t, n) =1 - a/ — du 

rni r(n+l) 

Also using the result 


(1.49) 


M e“^ 

1 - / dz 

0 r (n+l) 


(1.50) 


Putting M = (t - nn) on both sides of (1.50) 


^ V J A (t-n7r) ] 3 

2 e“ A(t-n7T) 

j=0 j ! 


^(t-nrr) ^-z 

-1 - J dz 

0 r (n+l) 


Further on substitution of z = A (u-nrr) 


z = 0 


u = nv 


A (t-HTl) 


u = t 




A ■ xt , 

1- J e~ A(^ [ A(u-n77)^] du 

f(n+l) n'/T 


(1-51) 
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COUNTER MODELS 


Comparing (1.51) with (1.49) 

" A ^-t- n^^ t A(t-n77) 

W (t; n) = S 

j =0 j ! 

P[X = n] = W(t; n) - W(t; n-1) 

j =0 j ! 

\ ri l\Tr\ A(t-(n-l)77) ]3 

S e“ A[l-{n-l) 7 r) J 

3=0 j ! 


(1.52) 


(1.53) 


and n<[t/77] where [t/'/r] refers to the greatest integer 
contained in t/n 


1.12 COUNTER MODEL TYPE I (RANDOM VARIABLE DEAD TIME):" 

We illustrate the concept of counter model 
type I with the following problem : 

Problem : A particle arrives at t=0 and locks 
the counter for a dead time of duration Y^. With the 
registration of the particle the counter is blocked for 
a time of length say ¥ 2 * The next particle to be 
registered is that of the first arrival once the 
counter is fixed. The process is repeated where the 
successive locking times denoted by Y^, Y 2 , Y 3 are 
assumed to be independent with a common distribution 
P[Yj,. < Y] = G(Y) 
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and independent of the arrival process. We obtain the 
waiting time distribution of the inter-»arrival of the 
process. 

Let Z = Y + Yy. 

Where Y is the dead time following the registration 
and Yy is the residual life time before the counter is 
locked. 

Where Y is the dead time and Yy be the residual 
life time before the counter is locked. 

Denoting by g(-) and f(.) the density 

functions of the dead time and the residual free time 
(before the system is again being locked) have 
P[Z<Y + Yy<Z + d2] =|5(z)d2 

2 

= J g(y) f(z-y) dy 

0 

1.13 COUNTER MODEL TYPE II (FIXED DEAD TIME): 

Here the locking mechanism, is more 
complicated. As before, an incoming signal is 
registered if any only if it arrives when the counter 
is free. In type I counter only recorded particles 
induced the counter to lock. For Type II Counter every 
arriving signal can prolong the dead time in the 
counter, the associated locking times being added 
concurrently 
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r 

. — * * * 

1 I 

* = Arrival (registration) 

I 1 = Dead time 

Figure : 

Let 7 -^ < 7 2 < registration 

times of the first, second and the nth renewal and 
arrivals take place with Poisson rate A- The length of 
the dead time following every registration is tt . 

We assume [Tj^ - independent . 

identically distributed random variables. 

P < X] = 1 - e” (1.55) 

P f^n ~'^n-l ^ X] = F'(x) (1.56) 

where the structure of F(x) is to be obtained. 

Let P [X < < X + dx] = f(x) dx (1.57) 

00 

0 (s)= /e"2t f(t)dt = L[f(t)] 

0 

Let u(t) = dU(t)/dt be the renewal density of the 
process and U (t) be the renewal function. 



* = Arrival (registration) 

I 1 = Dead time 

Figure ; 

Let < T 2 < T ^ are the registration 

times of the first, second and the nth renewal and 
arrivals take place with Poisson rate X- The length of 
the dead time following every registration is tt . 

We assume [Tj^ - are the independent . 

identically distributed random variables. 

P < X] = 1 - e" Ax (1.55) 

P [Tj^ “'^n-1 2 X] = F(x) (1.56) 

where the structure of F(x) is to be obtained. 

Let P [X < < X + dx] = f(x) dx (1.57) 

00 

<p (s)= /e“^^ f(t)dt = L[f(t)] 

0 

Let u(t) = dU(t)/dt be the renewal density of the 
process and U (t) be the renewal function. 


Then the ^L apla ce transform of the Renewal Density is 
given by ' — ~ — ' 


L[U(t) ] 


L[fi(t) ] 
1-L[f (t) ] 


00 

L[f^(t)] C e" dt 


A 


Also L[U(t) ] = L 


dU(t) 

■dt 


A +s 

__ X 
A + s 
1 - 0 (S) 


(1.58) 


1 - 0 (S) = 


A 


A + s 


00 


f -St 

3 [e d(t)] 

0 d(t) 


^ A 

A + s 


00 


J [e 

0 


-St 


dU(t) 

d(t) 


d(t) ] J ^ = 0(S) 


To obtain dU(t)/dt we proceed as follows : 

U(t+dt) - U(t) = e- At -f (5t) if t<77 

given that initially the counter is free. 


Also U(t+dt) 
dU(t) 
dt 


U(t) 


5 A<^t + 0(5t) if t>-n 


dU{t) 

L[ ] = 

dt 


A e At if ^ 

A e" A”' if f > 

J e- -t 

0 dt 


dt 


(1.59 


(1.60 


(1.61) 

(1.62) 

(1.63) 
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H- 


St At 


00 

dt + / e“ 


dt 


n 


AJ e A+s) - Avr °°r 

0 i 




L[U(t) ] 


A e~ A ( A+s) 


TT 


Ae (A'**s)7T ^ 

+ (1.63 ' ) 

A + S A + S 


U(t) = 


dU(t} 

dt 


+ 

A +s 


- i ^ 1 

A e A'ff ( j g-svr 

s A + S 


(1.64) 


-1 


A 

A+S 


A ® 


- Ay 


-1 


A e" Att^ 


1 1 


A+S 


= \ e A'n (2-e A tj 


,11 

A e~ An(__ ___^ g-sw 

s A + S 


= A e“ A'/t A(t-n-) j 

if t > 7r 


A e' 


A'n 


if t < rr 


9 by the se_cond_shi fting property of Laplace Transfor m. 

b(t) = A e + A e~ A'*' (i-e“ A(t--«)^ 

= A At + g- A'^ if z < n 


( 1 . 65 ) 



1 [dU(t)] 1 

L[U(t) ] = - L — L[U(t) ] 

s dt s 


A 11- 

/\ + s S ^ + S' 


.-( A-fcs)7r 


U(t) = (1-e At) + Ae" [(t-TT)- -{1-e" A (t-7r) ) 

A 

if t > TT 

= 1 - e~ At if t < n (1.66) 

1.14 COUNTER MODEL TYPE II (RANDOM VARIABLE DEAD TIME): 

This counter process is quite difficult to 
analyse in general form. However, the results are known 
only for process with Poisson inputs having arrival 
rate A - 

Let p(t) be the probability that the counter is 
free at a time t/i and registration is possible. 

It is to show that 


P(t) = e" A I [ 1 - G(y)] dy (1-S7) 

0 

where G(y) represents the cumulative density function 
of the dead time distribution. 
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Given n occurrences of a Poisson process in 
the interval (0,t] the distribution of occurrence 
times is the same as that of n independent random 
variables taken from a uniform distribution in (0,t]. 

Proof : The counter is free at time T = t if and 
only if all dead periods engendered by these signals 
have been terminated before t. 



(0,t] has a Poisson distribution with mean 


From the law of total probability 


00 1 t 

P(t) = S [- J G(t-y)dy]:i 
j = 0 t 0 


( ;,t) 3 e 


(At) 'J . e~ 

2%[— \G(t-y)dy]3 

j = 6 t h j! 


(1.70) 


00 t 

= e“ At s ^ f [G(t-y)dy]3 /j! 

j=0 0 


e At e A [G(t-y)dy] 
0 


A dy e A 


e A) [G(t-y)dy] 


= e“ A e AJ G(t)df 

0 


P(t) = e A / [l-G(y)dy] 
0 


(1.71) 


Counting v;ith our assumption / p(t) is the 


probability density. 


dU(t) 


/\ P(t) 


where U(t) represent the rene'wal function 


in (0,t] P (a signal appearing in [t, t + 6't)]. 



= P(t) <st + 0 (5t) 

U(t+(St) = U{t) + [1. ^ (5t)p(t)+0 (l-p /(St) ] + 0((St) 
liiTi U(t+5t) — U(t) 

t - 0 = ^P(t) 

t 

i-e. dU(t) 

- = 
dt 

U(t) = A P (T) dr [ U(0) = 0) ] (1.72) 

Putting (1.71) in (1.72) 

U(t) = A/ e- A [(1-G(y)]dy 

This proves the result. 
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CHAPTER 2 


APPLICATION OF PALM PROBABILITY FOR OBTAINING WAITING 
TIME DISTRIBUTION IN WEIGHTED POISSON PROCESS 

2.11 The distribution of = { Tj^ - n = 

1,2,3... where Tj^ represents the time for the n th 
renewal in a Poisson process with parameter ^ is given 
by the density function. 

f(t| A) = A e" ; 0 <t<oo; A>0 (2.1) 

However, if A» instead of being a constant assumes a 
probability distribution with density 

a^^ 

0( A) = e-^ A ;^k-l (2.2) 

r(k) 

then the unconditional distribution of time till the 
first renewal (or arrival) is given by 

uo k 

j f (tl A) )^ ( A) d A = rT7- (2.3) 

0 ' (a+t)^"^-^ 

• • • 1” h 

But the waiting time distribution for the n arrival 

given that the first took place at T=t is not 
certainly, the n-fold convolution of (2.3). Because of 
v;eighting, the renewal structure of the Poisson process 


is completely destroyed leading to interarrival * 
distributions having infinitely divisible structure but 
with dependent increments. Denoting by x(t) the number 
of events in (O.t] from a poisson process weighted by a 
Gamma distribution is given by (2.2) it follows that 
for- t > s Cov [X(s), X(t) - X(s)] 

= E [X(s), X(t) - X(s)| /] - [E( A)3^ s(t-s) 

= E [ s(t-s)] - (E(A)]^ s(t-s) 

= s(t-s) var ( A) ^ 0 (2.4) 

a^ 

Were 0(A) = e~^A ; 0 < A < 

r(k) a,k > 0 

00 

1 f 1 

E(- ) = J - 0(A) dA 
A " 0 A " 

00 

j 1 a^ e“®A dA 
0 A " r(k) 

oo 

a^ ^ e~^A dA 

r(k) 0 

r(k-2) 

r(k) 

= ; k>2 

(k-l)(k-2) 

f (tj_) <00 
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E(TiU) 


/ti e-At4 dti 

0 

00 

= A ^ -1 e-^t^ dt^ 

A r(2) 1 

A^ A 

This makes the problem of obtaining the 

probability distribution of the interarrival time in a 
compound Poisson process (or any such compound 
dependent process) somewhat complicated. On the other 
hand such interarrival time distributions are often 
considered as very useful for practical purposes; say 
while obtaining the distribution of interpregnancy (or 
interbirth) intervals or the v/aiting time distribution 
between two morbidity spells or the interarrival time 
distribution between two accidents. In all the cases 
the hazard rate A('t) (or the intensity) 'even for the 
same t varies between individuals. 

Cox and Isham (1980) have dealt with a series of 
such problems on dependent process classifying the 
problems according to the nature of dependence. This 
chapter is devoted to exhibiting an application of Palm 


probability to obtain the interarrival distribution 


between the first and the r^*^ arrival (r = 2.3....-) 
for a dependent process. 

2.1 THE RESULTS: 

Defining Palm probability as the conditional 

probability of r events (r = 0, 1, 2, 3 ) in 

(0,t] given than an event has occurred at time T = O, 
Khinchine (1960) gives 

k t 

= J [p^_^ (r) - dr (2.5) 

a 0 

and 

k t 

=1- — I pQ (t) dr (2.6) 

a 0 

where represents the unconditional probability of 

r events (arrivals) in (0,t] following a weighted 
Poisson process (weighted by Gamma ' distribution as in 
(2.2) and 

k lim [1 - 

t-0 (2.7) 

a t 

is the intensity of the process (Khinchine, (I960)] 


Then 


V, 


(t) = E (e- At/ ^ At A 


(a+t)' 


. . ( 2 . 8 ) 


Using (2.6) and (2.7) we have 
k 


ra 
-a+t-^ 


k t 

1 i 

a 0 ^ 


which on differentiation gives 

k+ 1 


ra 
-a-f-t-^ 


Po<^> 


(2.9) 


Again - P[T-j^ > t j an event has occurred at 

T = O] . Where is the random time of the first 

event following the occurrence of an event at T=0. 
Hence 

d (k+1) ■ 

fl (t|.) = -7., .. (2.10) 


dt (a+t)^^-^ 

where f 2 ^(t|.) is the conditional density of the 
v;aiting time of the second event given that the first 
has occurred at T = O. (2.10) is obviously different 
from the unconditional waiting time density of the 
first arrival given in (2.3). Proceeding in this way, 
putting r = 1 and differentiating both sides of (2.5) 
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we have 


1 


V't (t) 


= — 


d a t 

— _[K( — )k (— )] 

dt a+t a+t 


fi (t) 


a (a+t) 


,?Si(t)] 


(k+1) t 


(a+t) 


. (2 . 15) 


Then , 


Fp(t| .) 


i-j^o(t) - J^3^(t) 


( 2 . 12 ) 


where F 2 (t|.) is the conditional cumulative 
distribution function (c.d.f.) of the waiting time 
distribution of the second arrival given than at T=0 
the first arrival took place. Then, 


F2(t| .) 


(k+1) 


aJ^+2 (j,+2) ka^^+l 


(a+t) 


(a+t) 


k(k+2)t a 
(a+t) 


(a+t) 


ak+1 


( a+t ) 


ta*^"^^ (k+1) (k+2) 


(a+t) 


(2.13) 


which is the . interarrival density function between the 


first and the third arrival given that the first 
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arrival took place at T=0 in a compound Poisson 
process. Proceeding precisely in the same way we have. 


V>2(t) 


k 

a 

d p (k+1 ) a 
dt 


|52(t) 




[0i(t) - i 

S2(t)] 

(2.14) 

k t 1 

k r 

(K+l)a^'^^ 

J<( — )2 , 

a+t 

a - 

- Uo(t) 

(a+t) 

t2 gk+l 

(k+1) 

k+2 


k+3 


( 2 . 11 ) 


2! (a+t) 

1 - + f^{t) + f^it) ] 


1 - 


a(k+l) (K+i) (k+2) t^ gk+l t(t+k) a^'*'^ 


t- 2 ! (a+t) 
when we get 


k+3 


(a+t) 


k+2 


fp (t 1 . ) 


a^’^^ (k+2) (k+3)t2 


2 ! (a+t) 


k+4 


(a+t) 


k+1 


(2.16) 


using 


d 

r n-1 1 


1 - S (t) 

dt 

L ‘ J 

and the recurrence relation 

= fn-1 

(tl.) -jfin-2(t) + 
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(2.17) 


k 


n-1 


■ 4 ,- 


we have 


tn 1 gk+i (k+2) 

(n-1) ! 


(2.19) 


which provides the distribution of the time between the 
first and the (n+1)^*^ arrival given that the first 

arrival occurred at T=0 for a compound Poisson process 
weighted by a Gamma distribution. 

The treatment for other dependent processes is 
precisely the same as illustrated in the foregoing 
example . 

REMARKS : 

If (Tj^ - Tj^_^) is the v/aiting time for the i^^ 

order of conception (i = 2,3, ) or (T^ - Tj_ 2 ) tnay 

be called 'Interconception' interval and the conception 
rate ^ (even taken independent of i) varies from 

individual to individual conforming to some" probability 
distribution say Gamma distribution, then because of 

wa it ing ,\ . 

(i) Renev;al intervals (T^ - i = 2,3, 

v;ill cease to become i.i.d.r.v.'s (or i.d.r.v.'s) 

contrary to the traditional assumption. 


(ii) The num. ber of renewals even in two non 

overlapping intervals will be correlated. The extent of 
correlation because of the weighting of the process. 
The renewal structure is completely destroyed; leading 
to the process to conform to infinitely divisible 
distribution with dependent increments. 

CORRELATION BETWEEN T. AND ^ (j > i) : 

i J~i J ^ 

Let us assume the probability distribution of to be 

a^ 

0( A) = A A^"^; 0^ A 

r(k) a, k>0 

We have 

COV (Tj^, Tj)=E[E(Tj,Tj 1 A)]-E[{E(Til A)} {E(E(Tj | A)}] 
E(Tj_Tjl A) = E[Ti(Tj+Tj_il A)] 

E(Ti2| X) + [COV (T^, Tj_^| A) 

+ E(Ti| A)E(Tj_il A) ] 

Now COV (T^, Tj_^) = 0 

Since Tj^ and being tv.'o non overlapping intervals. 

i(i+l) 1 

E(Ti^l A) = 

2 

j-i 

and E(T._J /) = 

A 

ECl^Tj! A) = i(t>i) 


/ V X ^ j /X J 


A 


E(-;-; 
A- . 
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k>2 


since EC — -) = ; 

A (k-1) (k-2) 


i(j+l) a2 

ECT^T.) = 

(k-1) (k-2) 


i(j+k-l) 

Cov (T^T^) = ; k>2 

(k-l)2 (k-2) 


i(i+k-l) 

Var (T.) = ; k>2 

(k-1)^ (k-2) 

Cov (T^,T-j) = Cov (T^, + T-j_^) 

Next \-ie derive Cov (T , T j_ 2 ^) = Cov(Tj_, 'J-'i+T j _ 2 _) 

= E(Tj_, (Ti+Tj.;^) ) ■ E(Ti+Tj_2_) 

= E(TiM+E(TiTj_^)- [E(Ti]2- E(Ti 

= Var(T^) +Cov(T^ ,Tj_ 2 _) 

Cov(Ti,Tj_;^) = Cov(Ti,Tj) - Var(Ti) 

a^ i(j-l) 

= , k > 2 

(k-1) = (k-2) 

Also Var(Tj_ 2 ^) = Var(Tj) + Var(T^) - 2Cov(T^,Tj) 
a2 (j-1) (j-i+k-1) 

(k - 1)^ (k - 2) 


Thus Cor(T|,Tj_^) 


i(j-i) 

V(i(i+k-l)) V( j-l) ( j-i+k-1) 


v/hich is independent of a. 

Now assuming \ to follow 6 distn. 


We have (p{t) = 


6(a.k) (l+A) 


; o< A ^ 3'^'^ ^>0 


Thus Cov(Ti, Tj) = E[E(T^TjlA)] - 

E[E(TiiA)] E[E(T^|A)] 


Now we have to estimate 


ECT^Tj lA) ] = E[Ti(Tj+Tj_jA) ] 


= E(T2JA) + [Cov(T= i,Tj_JX) 

+ E(Ti lA) E(Tj_^lA) ] 

Where Cov (Tj^ , T j_^) = 0 

As Tj^ and Tj_^ being two non overlapping intervals and 


i(i+l) 


E(T2 •) 


and E(T._JA) 


E(Ti,T.j lA) = i(i + l) E ( 


Where E 


1 y 1 ' 

B(a,k) 0 A' + 


1 / ^ d 

B(a,k) 0 A" - (l< + 2) 


6(a-2,k+2) 


5(a,k) 

r(a-2) r(k+ 2 ) 


r(a) r(k) 

k (k+1) 

= . for all a>2 

(a-1) {a-2) 

i(i+l)k(k+l) 

Thus E(Tj,T^|A) = 

(a“l) (a~2) 


i(i+l) k(k+i) 

Thus Cov(T.,t^) E [ 

(a-1) (a-2) 

i ( i+1) k (k+i) 

therefore Cov(T;,T.-) = 

^ (a-1) (a-2) 


i j 

- ] E [ - 

A / 

1 

- ij [E( — 


where 




0(A) A 


1 


B (a , k) 


(1 A®' ^ 

0 X (i+A 


1 


15 ( a , k ) 


/ 


^(a-1) - 1 
(l+^)afk 


1 


B(a-1, k+1) 



1 r(a-i) k r(k) 


Cov(T£,T^) = 
Cov(T^/r..j_^) 

Where V(Tj_) 


=> V(Tj,) 


Thei'e tore 


(a-1) r(a-i) r(k) 
k 


a > 1 


(a-1) 

i(i+l) k{k+l) 


k2 


ij 


(a-1) (a-2) (a-1) 2 

= Cov(T^,Tj) - Var(T^) 

i ( i+1) k (k+1) ij k- 


(a-1) (a-2) (a-l)= 

E[E(T^2 3 _ E[E(Ti!A) ] 


Var{T4) 


i ( i + 1 ) 

E( ) 

A = 

i 2 + i 


i2 

E{ — ) 


A 


i2 E( — ) 

A^ 


i2 + i i2 k(k+i) 

E ( ) 

A 2 (a-1) (a-2) 

i(i+i) k(k+l) i2 k(k+l) 


(a-1) (a-2) 
ik (k+1) 
(a-1) (a-2) 
ik(k+l) 
(a-1) (a-2) 


(a-1 ) (a-2 ) 

I (i+l) -i] 


ov 


■■(T, ,T^_,} = Cov(Ti,T^) - Var(T3_) 


i(i+i)k(k+i) ij k2 ik(k+l) 

(a-l) (a-2) (a-l)2 (a-l)(a-2) 

k(k+i) ij k2 

(a-l) (a-2) (a-l) 2 


Cor (Ti,Tj_^) 


Cov(Ti, "^j-l) 
VVar(Tj^) VVar(T^_j^) 


Where 


Var(tj_^) = E[E(Tj_^2 l;^) ] _ E[E(Tj_^|^) 2] 


(j-1) j 

= E [ ] 


(j-l)2 
E [ ] 

A^ 


= j(j-l) E[ — ] - (j-l)2 E 

A== 


(j-1) 


[j(j-l) - (j-l)M E[ - ] 

A" 


Var(T- ^) 


= [j-1] E[ - ] 

(j-1) k(k+l) 
(a-l) (a-2) 


Cor(Ti,Tj_^) 


i2k(k+l) (a-l) - (a-2) ijk 2 
(a-l) V k(k+l) i ( j-1) 


In this case Conception rate A follov/s B 
distribution of second kind with parameters a and k 


then its expectation is given by 


5G 


that means average conception depends on the parameters 
a and k. Hence E(A,) is directly proportional to a 
whereas indirectly proportional to k. 

Our results shows that correlation between 
two non-overlapping intervals is higher when ^ follows 
ft Beta distribution as compared to the case when X 
follows Gamma distributions. In particular we can 
compare the correlations between and 2^*^ order of 

conception- i.e for i=l , j=2 when/ follows gamma 

1 

distribution Corr (T, ,T,-T t ) = — 

k 

and in case / follows ft distribution, 

[ (k+1) (a-1) - 2 (a-2)k3k 

distribution Corr (T., ,Tp-T, ) = — — 

(a-l)yktk+l) 



CHAPTER - 3 

COMPARISON OP COHORT FERTILITIES USING PALM 
PROBABILISTIC TECHNIQUE 

3.0 In this chapter a suitable methodology is 
evolved to predict and then to compare the future 
fertility performance. Given the time of occurrence of 
first birth as tj: and t 2 respectively of the two 
cohorts, whose age of marriage differs by a time lag t, 
it is proposed to analyse the level of f ecundability 
changes between two cohorts by the average number of 
births to these cohorts of women given that the first 
birth to the first cohort of women took place at t^ and 
the first birth of the second cohort of women took 
place at t 2 (t 2 > t^^) leaving the residual fertility 
span (T-t^) and (T-t 2 ) respectively for the two classes 
of women who complete their fertility span at the age 
T. 

The purpose is to see whether postponement of the 
age at marriage is effective in reducing the fertility 
status by using the palm probability technique. 

3.1 PALM PROBABILITY; 

Stochastic modelling involving prediction of the 
number of events during a fixed period of time (0,t) or 
the prediction of the waiting time distribution between 



two consecutive events or between any two events say 
the ith and the jth event can be solved by an entirely 
new technique known as "Palm Probability". It is the 
conditional probability of a specific number of events 
given that the event has occurred at the beginning of 
the interval. 

Let us consider two consecutive intervals of 
length s and (t-s) respectively as shown in the 
following figure : 


! 1 I I 

j * I I 

0 s S (t-s) t 

Figure 3.2 

Further let the interval (0,s) be so small such that at 
most one arrival of the event with compound Poisson 
inputs (i.e. Poisson input weighted by a two 

parameter family of Gamma distribution for taking the 
variation of fecundability between individuals) is 
possible and let us assume that ther e is only one 
arrival in (0,s) without any loss of generality. Then 
as s - O, the point of arrival denoted by X (in Figure 
3.2) in the limiting position will be located in the 
beginning . of the interval (s,t) . What .is the 
probability of having a finite number of arrivals of 
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events, say k number of arrivals (k=0, 1,2, . . . . ) in the 
interval (s,t) of any finite length given an arrival at 
the beginning of s. What is the waiting time 
distribution of next arrival or the nth arrival in 
sequence follovjing that event (i.e. an event has 
occurred at s)?.' These limiting conditional 
probabilities are known as Palm Probabilities. 

Let (f)-^{t) = The conditional probability of k 

number of births in [0,t] given that first birth has 
occurred at T=0 [0j^(t) is a plam probability measure V 
k=0 , 1 , 2 , . . . n] and 

Vj.{t) = The unconditional probability of k 

births in (0,t); k=0,l,2, n and let the births 

occur with Poisson intensity where ^ varies from 

individual to individual following Gamma distribution. 

Then Palm's integral equations connecting 
and 0]^(t) are as follows (Khintchine (1960) 

t 

Vj.. (t) = ^ 

0 

k = 0, 1, 2, n (3.1) 

and 

t 

VQ(t) = h(t) = h(t) 1 (p^{t)dt (3.2) 


6U 


where 


h(t) represents the intensity of the hazard 


rate of having a birth at the age t. 

Since the hazard rate h(t) is a increasing 
function of t for a given woman with f ecundabi lity 
level ,A write 

h(tj “ A where s > 0 

as the conditional hazard rate. 

It follows that s = d loqh(t) 

Jit 

Where s is the time elasticity of the hazard rate i.e 
the measure of the percentage increases in the hazard 
rate with 1% measure in time or age and h(t)/t = o = A- 

For given we have time dependent poisson 

[A 

process V^(t|A) = dt [ e dt] 

k! 

v;hich is obtained by solving the differential equation 
V^(t+A(t)) - Vj._^(t) [ A + 00(t)] 

+ Vj.(t) [ l-Ae^'*^ +OA(t)] 
being the probability of occurance of k 
conceptions or events in time (0,t]. 

V (tU) = Prob. [no birth upto time t] 



"A A(t) 


= e 


where A(t) == 1 
S 

^(t) = v^ctlA) «#>{A) dA 

00 

= /e-AA{t) gk g-aA ^k-1 

0 — — 

r(k) 

= . f e-a+^(t)^ ^k-1 

r ( k ) 0 


( a -t A ( t ) ) ^‘ 

1 

/here A(t} = — - 1) 

k 

1 - [ — ]*" 


a +A(t) 



(-k) (a)^ 


_k+l 

cl 

-k 


Since lim A(t) = lim [ — [ ] 

t-^O t-*0 <5 

= 0 . 

lim A' (t) = - = -1. 

t “^0 

Which is independent of t showing the stationarity 
of the process. 

Also 

d Vo(t) a 

= d[ 

dt a + A(t) 

-k a^ A' (t) 


(a+A(t) 


By differentiating Palm's integral equation (3.2) 

d k 

— y^it) = - — cp^it) 

dt a 

a ciVQ(t) 

=> 00 ^^) "" " — 

k dt 




+ ^ A' (t) 


(a+A(t))^+l 


6 3 




\a-Tl 

6 


r where A' (t) = - e^^ (3.8') 

Since (Vj.j (t) j^) conforms to a differential 
difference equation given by n = 0 , 1 , 2 n . 

We have 


V„(t+5t| A) = V^_^(t) [h(t| A) + o(5t)] 

+ V„(t) [l-[h(tl A) + o{5t)]] (3.10) 

d 

— V^(t| A) = A [v„_^(tl A) - Vn(t| A)] 

dt 


We have 

Vn(t| A) 



dt 


(3.11) 


Again by Palm's integral equation, 
for k = 1. 

— Vi(t| A) = A [Vo(t| A) - V^(t| A)] 
dt 


k + 1 

V^(t| A) = -- A(t) [ 

a a + A(t) 

Where A(t) is given in (3.6) 

a d 

0^(t) = 0o(t) V^(t) 

k dt 

(k+l) A(t.) A’ (t) 

1 ! [a + A(t) 


(3.12) 


(3.13) 
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Precisely in a similar way 


V2(t) 


k (k+l) A(t) 2 a^ 
2 ! [a + A(t) 


<p2 


(k+l) (k+2) A' (t) [A(t) 


2 ! [a + A(t) ]■ 


( 3 . 14 ) 


V3(t) 


k(k+l) (k+2) a^^ [A(t)]^ 
3 ! [a + A(t) 


( 3 . 15 ) 


03 (t) 


(k+l) (k+2) (k+3) A'(t) (A(t)]^ 


(3.16) 


3! [a + A(t) ] 


Proceeding in this way we get 


0n(t) 


(k+l) (k+2) (k+n) a’^^^ A'(t) [A(t)] 

[a + A(t) 


( 3 . 17 ) 


Expected number of births in the residual fertility 
span (T-t^^) given that the first birth has occurred at 
t=t^ for the first cohort is 


S n (T-t^) 

n=0 


A(T-ti) -|_(j^^2) 

a3_+A(T-t3) -I 1 (^1 + A(T-t3) 


= f^ [a^, k^, 5, (T-t^)] 


( 3 . 18 ) 


M-e- 


and the expected nuniber of births in the residual 
fertility span, given that the first birth has occurred 


at t=t 2 for the second cohort is 


n=0 


n (T-tj) 


A(T-t,) 

1 

a2+A(T-t2) 


-(K +2) 12 
X 


K+1 


= fp [a^, kp, 5, (T-t.) 


A(T-t2) (K 2 +I) A' (T-t2) 
(a2 + A(T-t2) 

(3.19) 


where (a^^, , i = 1, 2 are parameters of the Gamma 
distribution corresponding to the cohorts I and II 
respectively. 

Therefore for a change of age at marriage by 


= ^ 2-^1 


for 


t2>ti 


there is change in the average level of births 


E = 


^l^^l' ^1' ^1' ^ 

^2' '^“'^2' ^ 


( 3 . 20 ) 


f^ > f 2 indicates fall in the overall fertility level 
f^^ < f 2 indicates increase in the fertility level; and 
ff ~ ^2 indicates the stationarity in the fertility 

level . 

3 . 2 ESTIMATION OF THE PARAMETERS OF THE MODEL : 

VJe have for the first cohort E^[X(T-t^)| there is 
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f 

i- 


a birth at t=t^] 


ACT-tj^) 

1 

. a 3 ^+A(T-t^) 


ll_ 1-(K +2 

-t. ) 1 


A(T-t^) (Kj^+l) A'(T-t^) 

(a^^ + A(T-t^) 


and for the second cohort 


say 


(3.18) 


E 9 [X(T-t 2 )l there is a birth at t-t 2 ] 


A(T-t2) l.(K +2 


3 i,- ] 


(a 2 + A(T-'t 2 ) 


= f 2 [32, ^ 2 , 5, (T-t 2 )] say (3.19) 

Then the variance of the number of births for the first 

cohort 

Var[X(T-t 3 ^) 1 birth at t=t 3 ^ = 


0 


2 n^ (T-t^) 

n=0 


0 


2 lx 02 ^^^ (T-t]^) 

L n=0 


(k^+Da^^'^^ A'(T-ti) MT-t^) 
a^+A(T-t^) 


(a3^+A(T-t^)-l 


(k +2) 

X 


-f 


1 


+ 1 (T-t^) (K^^ + l) (^2^+2) r A(T t^) 


a +A (T 1 j 


L (a-^+A(T-t^) 


(a. +A(T-t^)- 


(k 4-3) 
Z- 
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2 


ACT-tj^) (kj^+l)A' (T-t^^) 
aj^+A(T-tj^) 


___!i |-(k.+ 2 ) 

(a^+A(T-t^) ‘ 


f. [a,, k,, 5 , (T-tp)] 


say ( 3 . 21 ) 


and the variance of' the number of births corresponding 
to the second cohort 


n =0 


<P 

S n 
- n =0 


(T-t 2 ) 


(k2+l)a2^2^^ A'(T-t 2 ) r a. 


a 2 +A(T-t 2 ) 


(a 2 +A (T-t 2 ) 


-(k+ 1 ) 


+ 32 ^ 2 ^^ A' (T-t 2 ) (K 2 +I) (k 2 + 2 ) 


A(T-t 2 ) 

- (ai+A(T-t2)J 


■ 

- (a 2 +A(T“t 2 )' 


-(k + 3 ) 
3 . 


a2K2+l A(T-t2) (k2+l) A' (T-t2) 

a 2 +A (T-t 2 ) 




(a 2 +A(T-t 2 ) 


(k + 2 ) 

Z 


= [ 32 , k 2 , 5 , (T-t 2 )] say ( 3 . 22 ) 

Now by the method of moments ( 3 . 18 ) and ( 3 . 21 ) for 
given <S and T-t^ will enable us to estimate a^ and' 


by successive approximations. . 


Thus by the same procedure (3.19) and (3.22) for 
given 6 and T-t 2 the estimates of a 2 and ^2 will be 
obtained . 

On substitution of the estimates a^, a 2 , k 2 

and 5 for given (T-t^) and (T-t 2 ) one can examine the 
increasing or decreasing trend of fertility as observed 
in two cohorts who differ by t years at age at 
marriage. The inherent difference in the fertility, 
apart from the availability of total marital exposure, 
if any, is expressed in the estimated values of the 
parameters (a^, k^^) i = 1, 2 for the two cohorts. 

3.3 A NUMERICAL ILLUSTRATION: 

When t^ = 16, t 2 = 18, T-^ = 45 
Thus 

E^[X(T-t^) 1 1^] = f^(a^,k;^,5,T-t;^) = (A)=3.64,4.48,5.50 
Var 3 _[X(T-t 3 _) I t 3 _]=f 3 (a ^ , k ^ , 5 , T-t ^ ) = ( B) =10 . 49 ; 9 . 6 5 , 10.69 
T-t^ = 45-16 = 29 years 

E2[X(T-t2) it2] = f2 (a2,k2,6',T-t2) = (C) =3 . 09 , 4 . 40 , 4 . 62 

and 

Var2 [X (T-t2) | t2 ]=f4 ( 32 , k2 , 5 , T-t2 ) = (D) =9 . 4 3 , 15.20,7.90 


45-18 


2 7„ years. 


Solving (A) and (B) for the given set of hypothetical 
values as given above by starting with an approximate 
initial solution of a^ = a^° and we get by 

successive iteration 

a^ =0.35, 0.40 = 1.0, 1.5, 2.0 

Similarly starting v;ith an approximate initial solution 
of a 2 = 92^ ^2 ~ ^2*^ solve (C) and (D) by 

successive iterations and get 

a2 = 0.35, 0.35 k2 = 1.0, 2.0, 2.0 

Substituting the estimates of k^ and a 2 , k 2 in 

(3.20) we get three values of 

E = 0.27, -0.24, 0.48 respectively 
which shows that in the hypothetical exercise the 
first cohort has the fertility status not significantly 
different than that of the second one. 

3.4 TREND OF THE PALM PROBABILITY DISTRIBUTION (NUMBER 
OF BIRTHS DURING FIXED MARITAL EXPOSURE) 

A comparison of the mean and the variance of the 
Palm Probability distribution given in (3.17) shov;s 
that 

Variance = Mean + Cj^ (Mean) - (mean)^ 


7 0 

I 

■ :f 


v/here? 


(1 


e 


a 


(k+3) 


[a + 


5(T-t )) 




k'f2 

( — ) 


-- (1- e 5(l-e-<S(T-t^ )) j 

5 ■ 


i = 1, 2 


If k, 6 and a satisfy a relation so that Cj^ = 6 (mean) . 

For the non-negativity of the variance it is 
necessary that © > 1 for Mean > 1. 

= Variance = Mean [1 + - (mean)] 

Variance = Mean [1 + {mean (0-1}] 

In particular, it may be noted that if mean is greater 
than one, if 0 < © < 1, Variance < Mean; or © = 1 
=> Variance = Mean. 

0 < © = Variance > Mean. For 0=2, Variance = Mean = 
(Mean) ^ . 

Thus in this flexible model by suitable choices of 
the parameters, the mean/ variance ratio can be made in 
the same order as may be observed in the sample taken. 


CHAPTER - 4 

THE WAITING TIME DISTRIBUTION BETWEEN THE FIRST AND 
SECOND ORDER CONCEPTIONS USING PALM PROBABILITY TECHNIQUE 


4.1 Biswas and ■ Pachal (1983) have extended Singh 
(1964) model on a probability distribution on the 
time of first birth in which it is assumed that a 
conception taking place during the infecundable period 
following the effective marriage with conception rate 
As it is further noted that the fecundability 
parameter X varies from the individual to individual 
even with same group. Biswas and Pachal (1988) have 
considered weighting of the Poission Process by Gamma 
distribution the renewal structure of the process is 
completely distroyed leading to inter-arrival 
distribution having infinitely divisible structure with 
dependent increments. In our study we have considered 
weighing of the Poisson process by weibull distribution 
and the results so obtained has been compared with the 


earlier existing results. On 
fecundability parameter /\ follows 
with probability density function 


k 

0(t) = — 

a 


A-M 

a 


] 


k-1 


e 




assuming that the 
Weibull distribution 
given by 

k 

a > 0 , k > 0 . 
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0 (A) = 


0 < A ^ 


Jf(tlA) 0(X) 

} 

fx dX 


A-M 

Put ( ) = u 


i dA = du 


; A'*' e du 


)]f U=0 , cl-- 1 

m 

= Ja ^ ^ ^ 

0 

= /" e-At dA 

$ 

y considering the distribution for / as MeibuU the 
alculation becomes »re complicated tHerefore 
:onsider the distribution A as a Standard «eibull 


istribution. Hen 


ce Weibull distribution- is not 


esirable for the analysis in our ^tudy 


Feller (1957) 


showed 


taht the number of 


conceptions 


to a couple during 


time 


interval (o,t] 



follows a Poisson Process with 


e At 

P[ N(t) = k] = 

k! 

Where N(t) is the number of conceptions during the time 
interval (0',t]. 

The probability of no conception is P[N(t)] = e^At 
or P[X < t] = F(t) = 1 - eAt fot given The 

conditional waiting time distribution for the first 
conception is thus given by 

f(tjA) = A e“At . 0 < t < 00 

The distribution is truncated at T' under this 
assumption of finite reproductive span for a woman is 
given by 


f (tlA) 


A® 


-At 


J A e“At 

D 


1 - e“At 


; 0 < t < if ' 


and the unconditional distribution of the waiting time 
from marriage to first conception with the fecundable 
parameter following Pearsonian type III is given by 


f(t) = 

0(A) = 


k a^ S 


^-o (a+t+k 


a e 


-A 


A > 0. 
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Hence the unconditional distribution of t is given by 



= / f (tlpv) cp{\} dA 

= J \ a e dA 

= ' dA 

= . a/A^-l e-(^+t)A dA 
a r(2) 

(a+t)2 

a 

= ^ — ; a > 0 

(a+t ) ^ 


However , there are several ways in which the 
model can be generalised. Let us take h(t|A) = A ® i 
6 > 0. instead of h(t|A) = A based on which (4.5) is 
developed. 

Let G(z,t) = Probability distribution function of 
Vj.(t) for all r = 1,2, ...,n. = S 

a. 

= 2 J v^(tlA) 0(A) '^A 

rs# * 

-t c 

dt) -A/e~*^^dt 

where V (t) = — — e 

r ! 


A" 


e 5 

<S 


[ A A(t) -AA(t) 

G 


r 1 
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where A(t) 


(1 


- 


^ ^ (AA(t))^ -AA(t) a e"^A 

Thus G(z,t) = S e 

r ! 


dA 


0- (A(t) 


r(r+l) 


r ! 


(a+A(t) ) 


r+1 


z"- (A(t) r(r+l) 
= z 

ri (a+A(t))^+l 

« z^ [A(t)]^ r! 

= Z , — — 

r! [a+A(t) 

1 - Z A(t) 


0 ^ 

2 


a+A(t) a+A(t) 


]• 


11 

a+ A(t) z A(t) 

[ 1 ] 

a + A(t) 

1 

a-l- A(t) - ZA(t) 

the choosing of the Pearsonian type III distribution 
for the fecundability parameter X is that for suitable 
parameters the distribution is unimodel and has the 
value zero at the limits of its range zero and 
infinity. Hence it is reasonable to assume that most of 
the women in apopulation v/ill have expectations near 
the average and that the number of eacli fertility level 


will decrease to a small proportion of the total as the 
expectation of the child bearing fails to zero or 
increases to very high rates and thus rates and thus 
justifying the choice of type III function. The choice 
of this function is justified by its flexibility too. 

The expected waiting time for the first conception 
is given by 

a a ^ 

E{t) =[ — [ 1+ _ ] ] 

k-1 a+t' 

Also we have 


a 

E(t2) = [ 

(k-1) {k-2) 


and 


a ka 

Var(t) 1 

(k-1) 2 (k-2) 


2a -(kt'-2a) ( )^ ^ 

a+t' 


- ( 


4 

Aft' 


) 


k-1 


kt ’ ( k"* 1 ) +2 a 

( 

(k-2) 


aM^-^)] 

a+t' 


The maximum time exposure for the first conception 
from effective marriage is taken to be 5 years since we 
are concerned with the proportion of women who has 
their first birth during first live years of their 
marriage . 'Effective marriage ' in means that the wife 
was formally marr'ied, had begun living with her husband 
and had begun mensturat ing . Under the assumption that 


reproductive span i.e 120 months the women who have not 
concieved within five years of her marriage will 
concieve in next five years of their marriage. 

j 


k ^ 
= a^ Z 


(a+ j+kt ' ) ^ 


] 


On substituting the value of a anf k will provide 
us the probability of the first conception for a 
married women at the above mentioned age group. 


4.2 PROBABILITY MODEL FOR WAITING TIME OF SECOND 
CONCEPTION WHEN FIRST TOOK PLACE (AT T=0).:’ 

If the hazard rate for conception is given by 
h(t) = A e-5t which represents the decreasing 
function for 5 > 0 and A being the intensity function 
for conception. 

For given X", we have for time dependent poisson 
process 


Vn(tlA) 


t f St 

[X / e'^^dt] ■ tA J e dt] 

? [ e 


Vj^(t) being the probability ol occurance of n 




78 


conceptions or events in time (0,t]. 

Thus conditional probability of concieving upto 
time < t, given ^ is as follows: 

1 - Vo(t|A) = 1 - e o dt 

-AA(t) 

= e 

1 

where A(t) = (l-e*^^) 

6 


/V varies from individual to individual following 
Pearsonian type III . 

The probability that one conception takes place is 

f’o(t) = ^ 3 

d.t 

f (t) = j— dt 

(a+A(t) 

k e^t ^ 

f(t) , 

(a+A(t) ) ) 5a 
1 - ( 

1 + (S a 

Thus the expected waiting time from marriage to 
first conception, is obtained from 


E{t) = /t f(t) dt 
* 

and t' represents the waiting time of the second 

conception given the first has occured at t=0. 

1 - cp'(t) is the conditional probability of not 

concieving given first has occured at t=0. 

Thus 1 - <p'(t)} = f'(t) which is conditional 

ett 

probability density function of waiting time of the 
second conception given the first has occured at time 
t=0 . 

d a^'^^ A' (t) 

= _[ 1 _ ] 

dt (a+A(t))^*^^ 

where A'(t) = e^^ 

and A(t) = — ( 1 - e*^^ ) 

<5 

Thus the expected waiting time for the second 
conception given that one has occured at t=0 is given 
by 

00 

E(t * I t=0) = / t f ® (t) dt, 

0 

Since the reproductive wastage is less and the 
medical facilities are good thus the actual foetal 
wastage is less leading to certain interesting 
methodology describing the probability distribution of 
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the waiting time from the marrige to first conception 
and from first conception to second conception with 
time dependent increasing hazard rate for conception . 
Thus the present study will reveal that the expected 
waiting time from the marriage to first conception is 
longer than that from first to second. 



CHAPTER - 5 


A STERLIZATION POLICY USING MULTISTATE MARKOV CHAIN MODEL 

5.1 The number of surviving children surviving is 
taken as the basis for Sterilization Policy for the 
reduction of the birth rate in a population. In general 
, a sterilization policy is implemented after a couple 
has a desired number of children on assuming a certain 
hypothetical survival rates operating in future. 
However, this method has no rationale to support : 
viz., the long term effect of mortality in deciding the 
number of surviving children may not be properly 
accounted for. An improved technique of assessing the 
long term impact of such sterilization policies is 
therefore to consider the joint effect of fertility and 
mortality which is the appropriate rationale for the 
present exercise. A Multistate Markov Chain- model has 
been developed corresponding to varying fertility and 
mortality intensities at different levels (States) of 
surviving children of the couple. The asymptotic 
probabilities of iiaving a fixed number of children have 


been derived. 


5.2 DEVELOPMENT OF THE MODEL : 

NOTATIONS; 

(i) ” Probability of n number of children upto 

a time t, (n = 0.1. 2. 3 ) 

[ (ii) An ® “ Probability of having (n + 1) 

I 

j children in (t + £ ) given that n number of 

I 

(Children are surviving at time t. (n = 

i 

0.1. 2. 3 ) 

(iii) + 0(ft) = Probability of having (n-1) 

children in (t + a ) given that n number of 

children are surviving at time t, (n = 

0 . 1 . 2. 3 ) 

(iv) M = A possible upper limit of the number of 

children to a couple. 

(v) (M-K+1) = Number of surviving children at the 

time of sterilization based on certain policy. 1< k 

. —t# 

<(M+1) 

P^Ct + O) ( 1-Ao^ +0(-5)) + ^ ..(5.1) 

Pj^(t+&) = Pn(t)[ l-(Xn+^n)^ 

+ Pn-i(t) [An-i^+0('^) ] 

+ Pp+l^^) f^’n+1^^ 0(a) ] .(5.2) 

for 0 < n < M-k 
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(5.1) => Po'(t) “ - >>0 PoW + «1 (S-'l) 

So on 

(5.2) => Pn* ^n-1 

+Pn+i(t) Mn+l 

for 0< n < M-k 


and 

(5.3) => P*M-k+l^^^ "" ■%-k+l ^ M-K+1^^^ 

■*■ /^M-k^M-k(^^ •'* 

5 . 3 STEADY STATE SOLUTION : 

steady state solution of (5.4), (5.5) and (5.6) 

are obtainable by assuming 

P (00) = Pj^(°°) = Pn-k+l^"^ ' 

Applying these conditions and solving 

successively the (M-K+2) equations, the parameters 

1 M Mo, •• parameters are 

Ai( •• AM-k+1' ^1' ''2' Mk 

obtained as follov^s . 

Ao 


Pn(°°) = -Hq = [ ^ 




XqH 

M 1 M 2 


XqAi- -Am- ( k+i) .,-1 
+ . . + — ^ ^ -1 
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Pn(°°) 


-^0^1 A 


^n-1 

V2 


M-k+l 


[ 1 + 


for all o<n<M-k 

%-k+l %-k+l %-k 

+ + . 

^M-k '-^M-k AM-(k+l) 


An-k ^M-(k+i) • • -Ao 


Examples : 


Case I : For M-K-1 = 2 


i.e. sterilization is performed on attainment of two 
surviving children. Parameter matrix is fXo ') 

^ M. M„ / 


Ml M2 


TTo = [1 + 


j-1 

II M 1 M 2 


TTo = [ 1 + 


M 2 M 1 


Ai -^0 


Case II : For M-K+1 = 3. 


i,e. sterilization is performed on attainment of three 
surviving children. Parameter matrix is / \q A 2 1 

M2 M3 / 


I 


^1 ~ Ao ^0 
^2 “ Ao Ai^^o 

7 T 2 = [ 1+ M3 +M3M2+M3M2Ml]~^ 

For presenting the effect of different levels of 
sterilization at various parity levels of the couples 
(from two to three surviving children), we define the 
fertility and mortality classes as follows: 

FERTILITY AND MORTALITY CLASS I ; 

Aq “- 20 , ~ .15, A 2 ~ • 

As ~* 10 » A 4 ~ .08, A 5 ~ -06; 

M 3 = 0.05, M2=M3=M4=M5=0.04 

FERTILITY AND MORTALITY CLASS II : 

Aq ~' 20 , ^3 = .25, A 2 ~ 

_A 3 =.20, A4 ~ .lo, As” .10, 

M 3 = 0.06, M2=M3=M4=M5=0. 05 
Steady state probability distributions of the 
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number 


levels 


of surviving children classified 
of sterlization is given below : 


by different 


Table 5.1 


Number of 
Surviving 
Children 

Levels 

2 

of Sterlization at 

3 4 

parties 

5 

0 

0.05 

0.02 

0.01 

0.002 

1 

0.20 

0.06 

0.02 

0.01 

2 

0.75 

0.23 

0.08 

0.04 

3 


0.69 

0.25 

0.11 

4 



0.63 

0.28 

5 




0.56 


Table 5.2 


Number of Levels of Sterlization at parties 


Surviving 


Children 2 


0 0.03 

1 0.16 

2 0.81 

3 

4 

5 


3 4 5 


0.01 

0.002 

0.001 

0.04 

-0.01 

0.002 

0.18 

0.04 

0.01 

0.78 

0. 19 

0.05 


0.76 

' 0.20 
0.73 



5.5 DISCUSSION : 


Proportion of surviving children at all levels of 
fertility and mortality classes under consideration 
decreases consistently when sterilization levels are 
relaxed. The proportion of couples having two surviving 
children in the long run is 75% when sterilization is 
made on attainment of two surviving children as against 
only 56% of the couples having five surviving children 
when sterilization level is relaxed to the extent of 
sterilizing couples only after having five surviving 
children. 

The same trend is more or less revealed in table 
corresponding to different sets of fertility and 
mortality levels as may immediately be seen on 
inspection of the diagonal elements of the matrix 
corresponding to the table. However, a better insight 
into the variation of the proportion of surviving 
children as a result of changes in the level of parity 
under consideration for adopting a sterilization policy 
may be obtained v/hen we consider the average no of 
surviving children under different F.M. set up. This is 
considered in the following table. 


Table 5.3 


Mean Number of Surviving children at different levels of 

Sterlization 


F.M. 

Class Level 

of Sterlization at number 

Children 

of Survivinc 


2 

3 4 

5 

I 

1.70 

2.60 3.49 

4 . 33 

II 

1.77 

2.73 3.69 

4 . 65 


Table 5.2 provides a clue for deciding the optimal 
level of sterilization depending on the mean num ber of 
surviving children. For example, if the policy is to 
ensure at least one serving male child in the long run 
(as it is most customary in India) per couple then it 
may be seen that irrespective of the variations in- the 
levels of the fertility and mortality the level of 
sterilization may optimally be made on the basis of 
three surviving children. Whereas if sterilization 
policy is relaxed to the extent of sterilizing couples 
only on the basis of more than three number of 
surviving children then one may expect a phenomenal 
increase in the growth rate ot the population not 
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consistent with the situation leading to adopt a 
sterilization policy motivated to cut down the growth 
rate of the population. On the other hand if 
sterilization policy is made on the basis of only two 
surviving children then the condition of ensuring at 
least one male child in the long run per couple may not 
be realized although the same may fulfil, the plan of 
reducing the growth rate in a fast increasing 
population. Such a policy may not be culturally and 
traditionally accepted by the people at large, thus the 
average number of children per family with one male 
child is three. 

5.6 MARTINGALE APPROACH : (STERILIZATION POLICY) 
INTRODUCTION: 

For a problem of sterilization of mothers with Iq 
number of suriviving children, attained by addition of 
a birth at the starting point of the observation, 
Assuming these mothers may be considered homogeneous 
with respect to age, i^ is a minimum number of 
children after v/hich the mother is eligible for 
sterl ization . 

However, if i^ is reduced to (i^ - 1) then the 
mother concerned, is discarded for sterilization once 
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* 1 , 


for all. Whereas those mothers whose numnber of 
surviving children move from i^ to (Iq + 1) first 

(without passing through any other state) is considered 
eligible for sterilization and they are sterilized at a 

random time S,j, on attainment of (i^ + 1) number of 

surviving children. Besides these two groups, the 
proportion of women who will remain in state number 

of surviving children throughout a period T from the 
starting point of the observation (the reference period 
7T which is considered large enough to cover all the 
three types of outcomes) will be discarded from the 
sterilization programme. 

The asymptotic estimator (for large T) of the 
above have been attempted to be obtained using 

Martingale stopping rule. 

5. 7 DEVELOPMENT OF THE MODEL: 

NOTATIONS: 

1) X(t) is the Number of surviving children at any time 
t of the birth and death process with birth and death 
parameters =^(i) = A i = M(i) ^ ^ respectively. 

.2) 7Tq = P[x(t) = Iq - l| X(0) = Iq ' 

X(j) = io 

for all j<t] 

where Wq is the asymptotic probability. 


3) a(T) - PlX(j) =iol*«» * io j£(0,T)). 

4) PK(t) ” P[X(t) = >^1 ; k > 0. 

We have 

P^(t+«) 

=Pv (t) ( 1- (Ak^+ +^ k-l (t)>^k-l^'^^k+l ^ 


for all k>0 

subject to P]<;_i(t) 

Ml 

Let us define 

Ai AiA2 


(5.8) 


+ . . + 


for all j>0 


M 2*^ . "M^ — 1 
^1^2**’* j-1 

. . . (5.9) 


corresponding to the parameters Ai A(i) ^i 

for* all 

Then it can Immediately be stated F(X(t)) is 

Martingale «ith respect to > 

T, = f(X(u, ; 0.uP.t (Karlin and Taylor (1P75„ . 

Let a stopping time be S^. 

= min(t It 0 X(t) = Xq ” ' -- 

where T=Srp .. (5.10) 

Where represents the stopping ^ 

• 4.1 nr X(t) =in”fO^ ‘ 

= in -1, or X(t) = lo+V°'' ^ 0 

0 • T is a Markov 

• . o - T and constant time T is a 

In this case S^, - i . 

riVicinrbinq states of tne 
time which are considered as ab 

j^(0) = in- optional- samp m 

process .given X(0) 0 ■ 


theor* 


f Martingale (Karlin and Taylor (1975)). 


E(f(X(ST))) = E{f(X(0))= E(f(iQ))= f(io), ...(5.11) 
Since X(0) =iQ, a fixed number by assumption. 

Then 

E(f(io)) = f(io) = E(f(XT)) 

= (I-ttq) [af {iQ) + (l-a) f (iQ+1) ] (f ( iQ-l) ) 

F{io) -af(io) 

= > ^ ^ 

f(iQ-l) -(l-a) f (I q+I) -af(iQ) 

(5.12) 

Then tTq proportion will not be sterilized and (l-a) 
(I-tTq) will be sterilized where 

a = e“ ( for large T (5.13) 

the reference period in the sterilization programme. 
PARTICULAR CASE: 

Let ig =2 

f(io-l) = 1 + 

Ml 

£'(io) = ^ ^ 

Ai 
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M1M2 




A1A2 



Hi 


[ 1^1 

^ o ^ t 


-t ’V)’')ll+Mi+!lL+i2 + 

A I 

[l+Hi+SSA] 

t 


A. 


‘ X. '2- 


»^4 />>) 

( 1+M2_+M;j_+M2) } 

A ^0. >' (5.14) 


Therefore, the total proportion 
be sterilized =[0!{l"'^o3 


)f mothers not to 


-(A - M )T foj. given T ... (5.15) 
where a = e tu y 


the expected stopping 


time Sn 


Next, to obtain 
i expe 

.lization, 100(1-.) 


^he expected time when policy decision 
which decides the exp 


taken for steri. 
while leaving 


ioo[a(i-^o3‘^'^o3 


mothers for 


sterilization programme. 


We have 


,(t, = ,(X(t))) 

- M(X(u))(G(X(u))-g(X(o)-l))<Su^^^^^^^_,^ 


is a Martingale over l°i 


random function g 
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where the expectation of Y(t) exists. Further, if g (i) 
be so chosen, so that 

\(i) [g(i+l)3-M(i)[g(i)-g{i-l)3 = 1 (5.17) 

for all i = 1,2,3.. 

then Y(t) = (g (X(t) ) -t) is a Martingale 

=> E(Y{S.p)) = E(Y(0))= E[g(X(t))3- E(S,pO ..(5.18) by 

Martingale optional sampling theorem. 

Further, Y(0) =g (X (0) ) =g ( ig) 

= q{io)= E(g{X(ST)))-E(SY)..(5.19) 

= E(g(X(ST)))-g(io) =E(S^). 

To obtain E (S.j,) , we require the solution of the 

difference equation given by (5.17). 

If g(i+l)-g(i) =^g(i)=V(i). 

We have the transfdormed equation of (5.17) as 

v(i) = (;,(i)V(i) - M(i)V(i)) = i. 

-‘ 4 , 

We may immediately get 

1 _ 

V(0) = 0,V(1) = (g(2) - g(i)) = 

Ai 

1 M 2 

V(2) = (g(3) - g(3)) = [1+ — ^ — ] 

X 2 <^1 

and in general 
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^{i) = 


[ 1+ 


i-± 




( 1+ _±- )(..) (1+ )] 

/\i-2 >^i-2 

(5.20) 


for i > O 

This gives the solution g(i) of (5.17) as 


1 1 
g(i)= — + 

Al -^2 


M 2 1 + n 

( 1+ — + ..••) ( 1+ -- ( 1+ — )) 

Ai ^^2 >1 


for i > 0 


and g(l) =0, 9(0) = 0 
In particular, for ig 2, g(io 


1 

g(3) = ■ 

Ai 

E(g(X(SY))) = 


A 2 

(l-Tig) ( 1 -a) 


7ro‘3(^o’^^^ * 


+a(l-vro)g(io)^ 


for ig = 3 

E(Srp) E(g(X(Srp))) - 9(io) 

Where f, g ana a are given in (5.9) 


(5.15) and 


(5.17) respectively. 


5.9 CONCLUSION: 

Thus by taking any 
be seen that for Iq " ^ 
may be the suitable most 


numerical Illustration it can 
(different values of M^, >i) 

. ■ ” .-unmiipr of surviving 

starting numnoet ut 


children in this Sterilization Programme under tte 
hypothetical values of the parameters. Thus this 
decides the policy of the sterilization programme. 
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CONCLUSION 

This Thesis pertains to the Analysis of some 
Stochastic Models and its applications in the area of 
fertility analysis and family planning programmes. 

The study is initiated with a technique based on 
Palm Probbability for the comparision of fertility 
statuses of two groups which are identical on Socio- 
economic parameters and different on age of effective 
marriage. The idea is to examine the policy of increase 
of conception by postponement of marriage. Thus it is 
believed that the methodology developed may be tested 
on a more exhaustive set of real data to derive more 
meaningful result. 

The model building exercise has been modified here 
by using Palm Probabilistic technique with restriction 
-of finite marital exposure. Thus comparision of results 
based upon the approach of infinite marital exposure 
with that of finite marital exposure ( Shreshta and 
Biswas 1985) reveals that the gain in precision in the 
estimate of interconception intervals for higher order 
of concptions. However for the higher order of 
conceptions further investigations are necessary in 
this case . 
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In chapter four the expressions for the expected 

waiting time for the live birth for all r and its 
variance were derived and a method to obtain the 
estimates of the parameters by the method of moments. 
Further the same method can be generalised by assuming 
p the probability of having a conception following a B 
distribution in confirmity with Potter and 
Parker { 1964 ) , which no doubt complicates the 
estimational techniques. But the complications remain 
worth adopting if we look into the practical utility of 
the same. Similarly fecundability parameter ^ may be 
assumed to vary among inidividuals following gamma 
distribution (with two parameters following BRASS 
(1958)) in which case the unconditional distribution of 
the waiting time of the first conception would be a 
Pareto distribution which will again complicate the 
estimation technique, but again complicatipn are seen 
in taking Weibull distribution. It has been noted that 
the estimation by the method of moments from Pareto 



distribution may be considerably simplified by . taking 
the correlation of interconceptive delays in a 
hetrogenious population. Another notable feature in the 


finding of the study is that when comparison is made 


between the estimates of 




the monthly probability of 
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conceptions using the probability models both with time 
independent hazard rate with those obtained with time 
dependent hazard, the probability values seem to be on 
the higher side in the latter case. A proportion of 
pregnancy wastage is less, and stress due to lack of sex 
education leads to under estimate of waiting time of 
conceptions and finally leading to the over estimate of 
monthly probabilities of conceptions. While developing 
the probability models for the waiting time 
distribution for different order of conceptions we have 
assumed that the parameters remain the same from parity 
to parity for the same individual i.e. for different 
order of conceptions the fecundability parameter 
follows the same distribution with the same pair of 
parameters a and k. This is leaving the situation in a 
dependent process where Palm probability technique has 
been made to obtain the waiting time distribution for 
the successive order of conceptions given that the 
previous order of conception has taken place at a given 
time T=0. Thus the effective marriage means that the 
wife was formally married and had begun living with her 
spouse. 

In Chapter-5 attempts have been made to develop 
Multisate Markov Chain model based on a density 
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dependent birth and death process to obtain the 
asymptotic probability ' of having fixed number of 
surviving children based upon the same distribution. 
Attempts have also been made to decide an optimal 
sterlization policy depending, upon the number of 
surviving children. If the policy is to ensure atleast 
one surviving male child in the long run then the 
results shows that the irrespective of the variations 
in the levels of fertility and mortality the levels of 
sterlization may be made on the basis of three 
surviving children. If the policy is relaxed to the 
extent of ster Using mother s only on the basis of less 
than three or more than three surviving children with 
atleast one male child , then one may expect the under 
or over growth rate which may not be consistent with 
the policy of reduction of growth rate. Finally, if an 
intermediate action based on the policy of three 
surviving children is adopted then the condition of 
ensuring one surviving male child in the long run per 
couple may only be realised with low probability. 
Martingales have been constructed on the continuous 
parameters ^ birth and death process to obtain the 
stopping time for sterlizing mother, for sterlization 
in the future, given that each of them have a fixed 


number of children. This methodology is useful in 
answering two basic questions 1- the cost of the 
sterlization programme and 2- stopping time. The 
percentage of mother's obtained for the sterlization 
programme based on the choice of different number of 
surviving children is indicative of the increase or 
decrease of the probable cost if we change over a more 
liberal policy while adopting the sterlization 
programme. Thus one can indirectly calculate the age of 
mother at the time of sterlization connected with the 
various sterlization programme. These two problems may 
be of great importance and should be kept in view in 
fram ing of the national policy programme. 

However, since the basic objective of this study 
is toestimate the correlation between waiting time for 
the first and second order of conception and its 
relation with the age of women and the application of 
Palm probability, Counter Model, Renewal theory and 
other Stochastic processes oriented methodologies as 
Martingale theory are some of the applications for 
arriving at the solution of fertility analysis which 
are methodological in nature and solutions are derived 
at theoritical plane and 'their justi f ication depends on 
more exhaustive real data. 
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